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Abstract 

We consider the totally asymmetric simple exclusion process 
(TASEP) with two-sided Bernoulli initial condition, i.e., with left den- 
sity p- and right density p+. We study the associated height function, 
whose discrete gradient is given by the particle occurrences. Macro- 
scopically one has a deterministic limit shape with a shock or a rar- 
efaction fan depending on the values of p± . We characterize the large 
time scaling limit of the multipoint fluctuations as a function of the 
densities p± and of the different macroscopic regions. Moreover, us- 
ing a slow decorrelation phenomena, the results are extended from 
fixed time to the whole space-time, except along the some directions 
(the characteristic solutions of the related Burgers equation) where 
the problem is still open. 

On the way to proving the results for TASEP, we obtain the limit 
processes for the fluctuations in a class of corner growth processes with 
external sources, of equivalently for the last passage time in a directed 
percolation model with two-sided boundary conditions. Additionally, 
we provide analogous results for eigenvalues of perturbed complex 
Wishart (sample covariance) matrices. 
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1 Introduction 



We consider the totally asymmetric simple exclusion process (TASEP) on 
Z. This is one of the basic one-dimensional interacting stochastic particle 
systems that, despite its simplicity, exhibits a number of interesting features. 
TASEP is a Markov process rj t with state space {0, 1} Z . For a given time 
t G M + and position x G Z, we say that site x is occupied at time t if rj t (x) = 1 
and it is empty if i] t (x) = (we can have at most one particle at each site: 
exclusion principle). The dynamics is defined as follows. Particles jump to 
the neighboring right site with rate 1 provided that the site is empty. Jumps 
are independent of each other and take place after an exponential waiting 
time with mean 1, which is counted from the time instant when the right 
neighbor site is empty (for a rigorous construction, see jl3|IHj). 

One may study a variety of slightly different observables of TASEP such 
as the total current, the location of a tagged particles or the TASEP height 
function. Here we focus on the height function, h t , defined from a TASEP 
configuration rj t as 



where N t is the total number of particles which jumped from site to site 1 
during the time interval [0,t]. 

In this paper we consider the simplest family of (random) initial condition, 
in which Shockwaves or rarefaction fans occur. More precisely, our initial con- 
dition is Bernoulli product measure with density p_ on Z_ = {. . . , —2, — 1} 
and p + on Z*j_ = {0, 1, . . .}. We refer to this as two-sided Bernoulli initial 
condition. Particular cases which have already been studied are: 

• the step- initial condition (p_ = 1 and p + = 0), where Z_ is completely 
filled. In this case, there is a rarefaction fan, the fluctuations of h t scale as 
t 1//3 , the correlation length as t 2//3 , and the limit process is the Airy2 process 
(see the case b = in p3]). The Airy2 process occurred first in closely related 
growth models [36, 49]. 

• Stationary initial condition (p = p_ = p + G (0,1)). The only station- 
ary and translation invariant measures are Bernoulli product measures with 
constant density p G [0, 1] (p = and p = 1 are however trivial) |42j . The 
scaling limit for the multi-point distribution of stationary TASEP has been 
recently unraveled in [BJ. 

Therefore we have only to focus on p + ^ p_ (the results below are the 




for j > 1, 
for j = 0, 
for j < — 1, 



(1.1) 
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content of Theorem 12.11 and are illustrated in Figure 12.11) . There are two 

cases: 

(a) p_ > p + . For large time t the asymptotic density decreases linearly from 
1 — p_ to 1 — p + over the region from (1 — 2pS)t to (1 — 2p + )t called a rar- 
efaction fan (see Figure |2~TT1 (a)). In this region the height fluctuations live 
on a t 1/3 and are governed by the Airy2 process like for step-initial condi- 
tion (with correlation length scaling as t 2 / 3 ). Around positions (1 — 2p±)t 
the randomness of the initial conditions start being relevant and there is a 
transition process from Airy 2 to Brownian Motion. When the fluctuations 
coming from the initial condition are on the t 1 / 2 scale, they dominate the 
fluctuations created by the dynamics (t 1 / 3 scale) and are governed by Brow- 
nian Motion. This is the case on the left and on the right of the rarefaction 
fan. 

(b) p_ < p + . For large time t there is a macroscopic shock with density jump 
from p_ to p + around the position (1 — p_ — p + )t (see Figure 12.11 (b)). For 
large time t, the fluctuations on the left and on the right of the shock are in- 
dependent. Of particular interest, is then the joint-distribution of the height 
function around the shockwave (1 — p_ — p + )t at different times. The initial 
conditions considered here are random and therefore looking far enough away 
the initial randomness becomes more important than the fluctuations created 
by the dynamics, which live on the t 1 ^ 3 scale only. For non-random initial 
conditions this does not happen and further limit processes arise, see [1^115] 
and, for one-sided random initial condition, see |16| . 

The results for the one-point distributions in Theorem 12.11 were conjectured 
in jH] and recently proven in [10]. The conjecture was based on univer- 
sality, since analogue results were available for a stochastic growth model 
(the polynuclear growth (PNG) model) [7], which is in the same universality 
class, named for Kardar-Parisi-Zhang (KPZ) [39]. The extension to multi- 
point distributions at fixed time in the PNG model was carried out in |33j 
(except for the case corresponding to stationary TASEP). 

Extensions away from fixed-time have been previously obtained in TASEP 
(with different type of initial conditions) in [T3J|33]. However, the extension 
was technically restricted to space-like paths, for which one could still get 
explicit expressions for the correlation functions. Our main result (Theo- 
rem 12. lj) is much more general, since it covers almost all space-time. In 
particular we can analyze situations where the correlation functions are not 
explicitly known! 

The only directions where the question of the limit process remains open 
are the characteristic solutions of the Burgers equation associated to TASEP, 
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also called characteristic lines [221EI]- Along these space-time lines the ap- 
propriate scaling limit is different, because the decorrelation occurs on a 
much longer time scale compared with the usual decorrelation length. This 
because second-class particles follow (on a macroscopic scale) exactly these 
trajectories. This phenomenon was first proven in a PNG model [31] and it is 
called slow decorrelation phenomenon (recently proven in greater generality 
in [20], see Proposition 12.51 below for TASEP). 

Methods in the proof of Theorem 12.11 

The proof of our main result, Theorem 12.11 employs a combination of many 
of the state-of-the-art methods in the study of TASEP fluctuations (exact 
determinantal correlation formulas, the connection to last passage percola- 
tion, coupling methods, and slow decorrelation). In outlining our proof we 
also provide a brief review of the literature on these different techniques. 

The first step in the proof is to establish a multipoint fluctuation result 
along a fixed space-time cut for a simpler initial condition corresponding to 
fixing p + = 0. For certain space-time cuts there exist exact determinantal 
expressions for the correlation functions. Specifically, if one considers the 
fixed-time cut for TASEP (i.e., the joint-distribution of the height at a fixed 
time) then there is a way to extend [TJJ[TB] to get the necessary formulas for 
the correlation functions [TT] . We consider a different cut which corresponds 
to a directed last passage percolation model with one-sided boundary condi- 
tion (see Section f2.2.2p . In that case the Schur process gives the multipoint 
correlation kernel, Equation ( 13. ip (for details on the Schur process and ap- 
plications see [TTlEZllinillZ]). In the TASEP setting, the Schur process is 
the process of a given (tagged) particle observed at different times. Using 
techniques of asymptotic analysis for Fredholm determinants we can extract 
our desired limit theorems from these formulas. This is done in Section 13.11 
and recorded as Proposition 12.41 As opposed to the related work of [33] , this 
is the only case for which we must appeal to the exact correlation formulas 
and take asymptotics. 

From this point on our proof relies entirely on probabilistic methods. 
The only case which is not covered by these other techniques is p_ = p + at 
the characteristic speed, but this was analyzed independently in [6]. These 
methods allow us to avoid the more involved shift and analytic continuation 
arguments of [B] in addition to the Schur process (see Remark 9 of [5J). 

The first probabilistic method we use is slow decorrelation (20] (given as 
Proposition 12. 51) . This implies that the fluctuation limit process for one-sided 
initial conditions extends away from cut on which it was proved (Propo- 
sition 12.81) . In order to bootstrap the one-sided process result to the full 
two-sided case we appeal to a coupling method introduced in [TO]. In fact, 
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versions of this method can be found in the literature in [31E21IM] under the 
names "microscopic Lax-Oleinik formula" or "strong monotonicity" . The 
new component offered by the coupling method of [10] is that due to a few 
(fairly simple) lemmas (see Section 13. 2j) one can now prove t 1//3 fluctuation 
results. The essential idea behind these methods is that statistics associated 
to complicated particle system or growth process initial conditions can often 
be written in terms of statistics associated to simpler particle systems which 
have been coupled to the original system. Once this connection is in place (for 
us Section 13. 3j) it is generally possible to translate asymptotic results about 
the simplier systems into results about the original (more complicated) sys- 
tem. This method is applied in Section 13.31 to prove Theorem 12.91 which is 
the last passage percolation equivalent of our main TASEP result. Finally, in 
Section 13.41 we show how to translate this back into a proof of Theorem 12.11 
One-side last passage percolation is closely related to the largest eigen- 
value of some complex Wishart (sample covariance) matrices |5j. Using the 
connection established in [HJ[22], we restate our one-sided last passage per- 
colation process result in terms of a random matrix eigenvalue process(see 
Theorem 12. lip . 

There are a variety of conjectured results which go under the title of uni- 
versality. The results of this paper deal with universality of the PNG and 
continuous time TASEP. However, TASEP is also the extreme case of the 
partially asymmetric version (PASEP), where particles can jump both left 
and right with different jump rates. For the one-point distribution function 
progress in this direction was made in in the early 1990s. Very re- 

cently, due to the efforts of Tracy and Widom [56-60J, Derrida and Gerschen- 
feld [21J, Balazs and Seppalainen [HE], Quastel and Valko [50J, Mountford 
and Guiol [15] significant progress has been made in answering this question 
in the general PASEP. Of particular note is the recent result of Tracy and 
Widom |59j which shows that the results of JTU] for TASEP with two-sided 
Bernoulli initial conditions extend to the PASEP setting for p + = and 
general values of p_. It seems hopeful that the integrable systems methods 
which proved useful in that paper will, eventually be able to deal with gen- 
eral two-sided Bernoulli initial conditions as well as multi-point distribution 
functions. With that eventuality in mind, this paper should serve as a guide 
in that pursuit. 

In addition to extending TASEP results to the context of the PASEP, 
Tracy and Widom's formula has played a prominent role in the long sought 
after calculation [HE3] of the one-point function for the KPZ stochastic PDE 
with narrow wedge initial condition (for an alternative approach using the 
replica trick see also [T9|l23]). 
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2 Results 

Here we present the limit results, first for TASEP, then for last passage 
percolation and we end with random matrices. The limit processes in the 
following statements are defined in Section 12.41 



2.1 Continuous time TASEP 

We want to analyze the fluctuations of the height function (11.11) with respect 
to the macroscopic behavior. Thus, the first quantity we need to determine 
is the limit shape 

/i ma (0:= Um~h t (m) (2.1) 

t->oo t 

which can be obtained by integrating the asymptotic macroscopic density of 
particles, f?(£,r), given heuristically by 

£>(£,t) := lim P(there is a particle at [£T] at time tT). (2.2) 

T— i>oo 

The average current of particles for a density g is g(l — g), thus g satisfies 
Burgers equation [5T] 

d T g + dt(g(l-g)) = 0. (2.3) 

The initial condition £>(£, 0) = p_ for £ < and g(£, 0) = p + for £ > gives: 
(a) for p-> p+, 

{p_ for£<l-2p_, 
(1-0/2 for£e[l-2p_,l-2p + ], (2.4) 
p + for i > 1 - 2p + , 
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(b) while for p_ < p + , 



(p + for £ > 1 - (p_ +p+). 

The characteristic lines^, (t,x(t))t>o, of the Burgers equation with constant 
density p are straight lines with speed 1 — 2p: {x(t)—x(0) = (1 — 2p)t, t > 0}. 
In the case of non-constant density (see case (a)), then all the rays leaving 
from the origin with speed £ G [1 — 2p_, 1 — 2p + ] are also characteristic lines 
(see Figure I2.2p . 

Translated into the limit shape using (11.11) . one obtains: 

(a) for p-> p+, 

'2p_(l-p_) + (l-2p_)£ for£<l-2p_, 
KJK) = { (l + a/2 for£e [l-2p_,l-2p + ], (2.6) 

2p + (l-p + ) + (l-2p + )£ for£>l-2p + , 

(b) while for p_ < p + , 

f2p_(l-p_) + (l-2p_)£ for£<l-(p_ + p + ), , , 
malU \2p + (l-p + ) + (l-2p + )£ for£>l-(p_ + p + ). 

This is illustrated in Figure 12.11 

For simplicity, we discuss the fixed time fluctuation results of our main 
Theorem 12 .11 which is however more general and holds for unequal times too. 
Concerning the fluctuations for p_ ^ p + , if we focus around a macroscopic 
position £t we have: 

Case (a) and 1 — 2p_ < £ < 1 — 2p + : the limit shape is curved and the 
behavior is like the one of step-initial condition, namely, for large time t the 
fluctuations scale as t 1 ^ 3 , correlations as t 2//3 , and the multi-point statistics are 
governed by the Airy 2 process, Ai- More precisely, there are two coefficients 
n h = (2(1 — £) 2 ) 1//3 an d «« = — (1 — £ 2 ) 2 / 3 /2 1 / 3 depending only on £ (compare 
with (12.91) below) such that for large time t, 

h(£t + TK h t 2/3 ) ~ £/i ma (£ + TK h r 1/3 ) + K v A 2 (r)t 1/3 . (2.8) 

Case (a) and £ = 1 — 2p_ (or £ = 1 — 2p + ): the influence of the randomness 
in the initial condition and the randomness built up by the dynamics are of 
the same order. The fluctuations of h t are on the t 1 ^ 3 scale with correlation 



1 These characteristics are the ones coming from the entropy condition [251161] . 
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<?(£,!) 



*?(£,!) 



-o. 



p- 




Figure 2.1: The asymptotic density g and the limit shape in the cases (a) 
p_ > p + and (b) p_ < p + . Transitions happen at £± = 1— 2p ± and shockwave 
at £ s = 1 — (p_ + p+). The fluctuations processes are also indicated: B and B' 
two independent standard Brownian motions, A% the Airy 2 process, */4bm_;.2 is 
the transition process from Brownian behavior to Airy 2 process, and ^2^bm 
is its time-reversed version. See Section 12.41 for definitions. 



scale t 2//3 and are governed by a transition process v4bm->2 between Brownian 
Motion behavior and the Airy 2 process. 

Case (a) and £ < 1 — 2p_ (or £ > 1 — 2p + ) or Case (b) away from the shock 
position: the influence of the initial randomness dominates and one has sim- 
ply Brownian Motion (with fluctuations scale t 1 ^ 2 and correlation scale t). 
The two sides are asymptotically independent. 

Case (b) at the shock position: the statistics of the height function is influ- 
enced by both the right and left particle densities. 

The stationary case, p_ = p + = p was analyzed already in [6] (see Theo- 
rem 1.7 therein). At £ = 1 — 2p the fluctuations are of order t 1//3 and there 
is a transition process *4 s tat over a distance of order t 2//3 to the Gaussian 
behavior. 
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In a related model [3T] the following slow decorrelation phenomenon was 
noticed: along the characteristic lines the height-height correlations live on 
a longer space-time scale than the fixed-time correlation scale. For instance, 
in the rarefaction fan, the height function at two space-time points on the 
same characteristic line, the first at time T and the second at time T + T U , 
v < 1, will differ by a deterministic factor (speed of growth x T u ) plus 
o^T 1 ' 3 ). This means that the two height functions (centered and rescaled by 
T -1 / 3 ) are asymptotically the same random variable (i.e., they are perfectly 
correlated on the T 1 / 3 scale). The proof in [31] uses several results of other 
papers and a considerable amount of work is needed to reproduce them for 
other models. While looking for a proof for the TASEP, we discovered a 
much simpler proof, which applies not only to TASEP but to a large number 
of models in the KPZ class, see [2D]- The statement for TASEP is reported 




in Proposition 12.51 

This allows us to extend the fixed-time statement to space-time. This is 
the reason for the following limit theorem: for £ E [1 — 2p_, 1 — 2p + ], let us 
set 



X(t, 9) = [£(T + 9T V ) + r(2(l - £ 2 ))V3 T 2/3j ; 
H{t, 6, s) = — ^— (T + ffT) + er(2(l - £ 2 )) l/ *T 2/3 + (r 2 - s) [i J f J T 1 / 3 



The value of H(r, 9, 0) is a generalization of the term th ma (^ + TKht 1 ^ 3 ) in 
(12.81) . namely the macroscopic approximation. Indeed, 



compare with ()2.6|) . while H(r,9,s) — H(r,9,0) measures the fluctuations. 
The definitions of the limit processes occurring in the following theorem are 
collected in Section 12.41 

Theorem 2.1. (a) Fix m E N, u E [0, 1), f E R, and p + E (0, 1], p_ E [0, 1). 

Then, for any choice of real numbers t\ < r 2 < . . . , r m , 0\, . . . , 9 m , and 
Si, . . . , s m , we have: 

(al) If p+ < P- and £ E (1 — 2p_, 1 — 2p + ), then 




(2.9) 




(2.10) 



lim P P\{h 

\fe=i 




T+6 k T»(X(T k ,9 k )) > H(T k ,9 k ,Sk)} 




P (f){A 2 (T k )<S k } 



(2.11) 
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(a2) If p + < P- and £ = 1 — 2p_, then 



lirn^P f]{h T+6kT ,(X(T k ,6 k )) > H(r kl 9 kl s k )} J = P f|M B M^ 2 (r fc ) < s fe } ) . 

(2.12) 

fa3j If p + = p- = p and £ = 1 — 2p, i/ien 



T lim P p){/i T+0fcT ,(X(r fc ,^)) > H(r k ,9 k ,s k )} = P f|{Atat(r fc ) < s fc + r fc 2 } 

(2.13) 

(b) Fix m G N, and p + G (0,1], p_ G [0,1). Fix rni,m s ,m r G and set 
m = mi + m s + m r and rea/ numbers 9\, . . . ,9 m . Consider a set of m space- 
time points with macroscopic coordinates (£j#jT, 9{T) . Let mi of the points be 
such that £j < 1 — p_ — max{p + , p_}, m b stzc/i i/iai & > 1 — p + — min{p + , p_}, 
and ra s on £/ie shockwave fa — 1 — p + — p_ i/ p + > p_/ T/ien 



Km P f){ ho kT (£ k 6 k T) > h m ^ k )9 k T - 2s k T 1 ' 2 } 
\fc=i , 

f| {B(e*(l-2p_-e*)(p-(l-P-)))<«*}) (2-14) 
fc=i 

(m 
n {i3 / (^(6+2 P+ -i)(p + (i-p + ))< Sfc } 
fc=TOi+l 

where B and B' are two independent copies of Brownian Motion. 

Remark 2.2. Although in the statement we fix 9i,...,9 m , the same 
holds true if they depend on T provided that they are uniformly bounded 
in T. What we need is that there exists a v < 1 such that 
lim T ^xJn(|0 fe T"|)/ln(T) < 1. For instance, we can take 9 k T v = 9 k T 2 / 3 
with 9 k fixed real numbers. 

Remark 2.3. The case p + < p_ and £ = 1 — 2p + can be recovered from (a2) 
by particle-hole symmetry. The entries in the Brownian motions in (I2.14p 
are (proportional to) the projections of the space-time points to time t = 
along the characteristics to the initial conditions; the proportionality takes 
just into account the variance of the random walk of the initial condition. 
This is illustrated in Figure 12.21 



The proof of Theorem 12.11 is in Section 13.41 It is a consequence of the 
corresponding result for last passage percolation (see Theorem I2.Q|) . together 
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t 1 1 




Figure 2.2: Illustration of the characteristic lines for (a) p_ > p + and (b) 
P- < p+. The fluctuations of the black points depend on fluctuation of their 
projections to the t = line. Points on the shockwave, the gray points, 
depend on the projections on the two directions. Finally, the fluctuations in 
the rarefaction fan (white points) do not depend on the initial randomness. 



with the slow decorrelation phenomena (see Proposition 12. 5p . The stationary 
case, (a3), was analyzed in [6] (see Theorem 1.7 therein); the r| term in (12 . 1 3[) 
compensates the fact that the scaling (12. 9p is not following the (straight) limit 
shape approximation. For presentation simplicity in [6J only the fixed-time 
result was stated, but slow decorrelation allow immediately to extend it as 
in Theorem 12.11 

2.2 Directed percolation 

In Section 12.2.11 we explain the precise connection between TASEP and last 
passage percolation (LPP). In Proposition 12.81 we give the asymptotics for 
one-sided LPP, which uses the determinantal structure of the Schur process 
and slow decorrelation. The extension to two-sided boundary conditions via 
coupling arguments is stated in Theorem 12.91 

2.2.1 Connection with TASEP 

We define a directed last passage percolation model by assigning random 
waiting times u>jj to each site in (Z+) 2 (Z+ = {0, 1, . . .jl. We require 
WijS are independent and exponentially distributed variables^ (to be speci- 
fied below). To every directed (up/right only) path tt from (0, 0) to (x, y) we 

2 We use the notation Exp (to) for a random variable which is exponentially distributed 
with mean to. 
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associate the waiting time T(ir) = Yl(ij)eir Wi >r Then, the last passage time 
from (0, 0) to (x, y) is the longest waiting time over all directed paths: 

L(x,y) = max T(tt). (2.15) 

It is well known that the height function for TASEP with our initial condi- 
tion is expressible in terms of a LPP model. Let us shortly recall this connec- 
tion which is established in full generality in [48j and briefly reexplained [6], 
extending the step-initial condition case considered in [33]. We label parti- 
cles from right to left and denote by x*.(i) the position of particle k at time 
t. We set the label so that • • • < x 2 (0) < xi(0) < < x (0) < x_i(0) < • • • . 

(a) For i,j > 1, Wij is the waiting time that particle j jump from site i—j — 1 
to site i—j (of course, the waiting time counted from the instant where site 
i—j is empty). Thus, Wij ~ Exp(l) random variables, i,j > 1. 

(b) On the other hand, the effect on the dynamics on x (t) due to the par- 
ticles on the its right is equivalent to set the jump rate of particle to be 
1 — p+ instead of 1. This is a consequence of Burke's Theorem [18]. Therefore 
we set to^o ~ Exp(l/(1 — p + )) for i > x (0) and w^i = otherwise. 

(c) By looking at the particle-hole transformation, we set wij ~ Exp(l/p_) 
for j > — xi(0) and w j = otherwise. Finally, we set to ,Q = 0. 

With this settings, the correspondence between last passage time, particle 
positions, and height function is the following: for Xk, Vk > 1, tk > 0, we have 

P {C% =1 {L{x k ,y k ) < t k }) = P (nT =1 {xyM > x k - Vk}) (2 lg) 

= P (p^ =l {h tk (x k - y k ) >x k + y k }) ■ 

Since x (0) ~ Geom(l — p + ) and —(1 + xi(0)) ~ Geom(p_), we can set 
Wi } i ~ Exp(l/(1 — p + )) for all i > 1, and wij ~ Exp(l/p_) for all j > 1 
without changing the large time asymptotics (see e.g. Proposition 2.2 in |32j) 
(but keeping w 0) o = 0). 

2.2.2 One-sided LPP 

As briefly mentioned in the introduction, the proof of our result uses a mix- 
ture of analytic and probabilistic methods. On the analytic side, we have to 
analyze the following LPP model (referred as one-sided LPP): 

I Exp(l/r7) fori = 0,j>l, (2.17) 
[o fori>0,j = 0, 



where rj G (0, 1] is a constant. We denote by L\ the last passage time for the 
waiting times (I2.17p . where 1 stands for one-sided. This problem is related 
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to TASEP with p + = and p_ = r\. Moreover, see Section 12.3} the statistics 
of Li are related to those of the largest eigenvalue of a perturbed Wishart 
(sample covariance) matrix. 

To set the scaling variables, we need an expression for the limit shape. 
Let us focus along the line y = 7 2 x. There are two cases: 

(a) for 7(1 + 7) -1 < 77 < 1, 



(b) for < 77 < 7(1 + 7)- 1 



lim L 1 (eT,e 7 2 T)=e(l + 7) : 



lim —LxiT, 7 2 T) 

T— >oo 1 



77 7] 



(2.18) 



(2.19) 



In the regime where the limit shape is (I2.18p . the fluctuations of L\ are 
of random matrix type with correlation on the T 2 / 3 scale and fluctuations in 
the T 1 / 3 scale. Therefore we introduce the scaling 



y( r ) 



1 



(1 + 7) 2 

,2 



T + 



2r 



(1 +7 )2/3 7 2/3 



^2/3 



7 



£(r, s) = T + 



I + 7) 2 
2r(l 



T 



(2.20) 



7) 



1/3 



2^2/3 



[S — T 



(l + 7 )2/3 



T l/3 



^2/3 ^1/3 

where the parameter s is a measure of the fluctuations with respect to £(t, 0), 
that is what we expect to see from (12.181) . Under this scaling, the height 
fluctuations are governed by the Airy 2 process, A2, up to the critical value 
r] = 7(1 + 7) _1 where there is a transition process, *4.bm^2, to the Brownian 
motion behavior. In the regime where the limit shape is (12.191) . the fluctu- 
ations will be governed by the boundary sources. They have fluctuation on 
the T 1//2 scale, correlation length of order T and limit process the Brownian 
Motion, B. This is precisely stated in following theorem. 

Proposition 2.4. Let A2, Am^2 and B be the processes defined in Sec- 



tion 2Jf_ 



(a) Fix m e N, ?) 6 (0, 1] and 7 G (0, 00) with 77 > 7(1 + 7) _1 . Then, for any 
given T\ < r 2 < . . . < r m and Si, . . . ,s m el, we have: 
(al) if i] > 7(1 + 7)" 1 , then 

/ m \ / m 

XimV ( p|{L!(x(7>),l/(7]fe)) < £(r k ,s k )} ) = P ( fl^ 2 ^) ^ Sk ^ 



vfe=l 



(2.21) 
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(a2) while ifr] — 7(1 + 7) 1 , then 



Km P ( ^[{Lxixir^yiTk)) < £{r k ,s k )} ) =P ( f)M B M^ 2 (r fc ) < s fc } 



vfc=l 



vfe=l 



(2.22) 

(b) Fix m G N ; 77 G (0, 1]. T/ien, /or any given 71 < 72 < • • • < 7m suc/i i/iai 
77 < 7i(l + 7i)~\ and si, . . . , s m G 1, we /iai>e 



& p (n{ £ iCr,Tjn< 




(2.23) 



This theorem is proved in Section I3TT1 using the Schur process and applying 
methods of asymptotic analysis. Having established this theorem for the one- 
sided boundary condition model above we use coupling methods to prove a 
general two-sided boundary condition theorem. 

To understand intuitively the cutoff 77 = 7(1 + 7) _1 there is a simple 
argument. The last passage path (the random directed path which achieves 
the last passage time) goes along the left boundary for some distance and 
then will depart. Restricting the set of paths to only those which go a certain 
macroscopic distance along the boundary and then depart into the bulk, one 
may use independence of the boundary and the bulk to establish a law of 
large number and fluctuation theorem for the restricted last passage time. 
If the mean I/77 of the boundary waiting times is large enough (77 small 
enough), the restricted law of large numbers will be maximized for a positive 
macroscopic distance along the boundary. In this case, the fluctuations will 
come entirely from the boundary fluctuations and they will be given by the 
standard CLT: Gaussian fluctuations on the scale T 1 / 2 . On the other hand, if 
the boundary waiting times are too small, the restricted law of large numbers 
will be maximized for a distance along the boundary of o(T 2 / 3 ), and hence 
the fluctuations will come from the bulk, which are known to be T 1//3 and 
GUE Tracy- Widom distributed. At the cutoff, the two fluctuations compete 
and yield a perturbation of the bulk fluctuations. This intuition will be useful 
in some of the arguments used in this paper. 

To extend Proposition 12.41 to points which have different y coordinates we 
use the fact that in certain directions (the characteristics) the last passage 
time fluctuations decorrelate not in the scale T 2 / 3 but rather in the scale of 
order T. This means that the passage time at two points at distance o(T) but 
on the same characteristic have the same fluctuations (up to o(T 1 / 3 )). This 



14 



phenomena, known as slow decorrelatiorjfl, was observed (and proven) in the 
related PNG model in [31] and then extended to a much greater generality 
within models in the KPZ universality class (20]. We recall the result needed 
here below. 

The characteristic lines for TASEP with particle density p move with 
speed 1 — 2p. In the present LPP picture this implies the following: 

(a) if 7 < yz - ; the line y = 7 2 x is a characteristic related to the rarefaction 
fan of TASEP, 

(b) while in the case 7 > yz - , the characteristic passing by (T, 7 2 T) is given 
by 

y = T 1 2 + (x- T)n 2 (l - T])- 2 . (2.24) 

With these preliminaries, we can state the slow decorrelation theorem for 
LPP. 

Proposition 2.5 (Corollary of Theorem 2.2 of [20]). 

(a) For 7 < define 



T 



1 + 7 ) 2 J ' L(l + 7) 2 _ 



7 2 T 



Q 



T + r 7 2 (T + r) 

t+th'L (1+7) 2 



Then, for any r ~ T u with v G [0, 1) and any given M > 0, it holds 
\im P (|Li(Q) - Li(P) - r\ > MT 1/3 ) = 0. 

(b) For 7 > define 

P=([T\, L 7 2 TJ ), Q = ( [T + r j , [ 7 2 T + rr/ 2 (l - r/)- 2 J ) . 
Then, for any r ~ T u with v G [0, 3/2) and any §roen M > 0, it holds 

> MT 1/2 



lim P 

T->oo 



L X {Q)-L X {P) 



0. 



(2.25) 
(2.26) 

(2.27) 

(2.28) 



Remark 2.6. In Proposition 12.51 7 can be also chosen to be T-dependent, 
provided that it converges to a fixed number in the T — > 00 limit. 



3 For flat interfaces as considered in [3S] the dynamic scale invariance [H] implies a 
scaling form for the temporal autocorrelation, from which one could expect to see a slow- 
decorrelation type of phenomenon. This phenomenon is also related with what is known 
as persistence, see [38] for KPZ class and [40 for Gaussian-type models. 
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Remark 2.7. We will use a generalization of Proposition 12.51 to joint dis- 
tributions. We let r be a vector r = (r 1; r 2 , . . . , r m ) with each ~ T Vi . 
Let us interprete P and L%{P) as vectors P = (Pi, P 2 , . . . , P m ) and 
L X (P) = {Li{Pi),Li(P 2 ),...,L 1 {P m )), and similarly for Q and Li{Q). Then, 
the theorem still holds, with absolute values replaced by Euclidean norms. 
Case (a) holds for max(^) < 1 and case (b) for max(z/j) < 3/2). This follows 
directly from triangular inequality the union probability bound. 

As immediate application of Proposition 12.51 (and Remark 12. 7p is the 
extension of Proposition 12.41 away from the fixed-?/ line. Indeed, often one 
considers the cut x + y = t and t is then interpreted as the time parameter 
in a stochastic growth model (see e.g. For that reason we consider 

the following modification of the scaling f)2.20p : for a given v G [0, 1), 

2 7 4/3 



X(T, 



y(r,e) 



i 

(T+7? 

7 2 
(l + 7 )2 

6T U + t 



(T + 6T U ) + r 



(T + 6T V ) - t- 



(1+ 7 2)(1+ 7 )2/3 

2 7 4 / 3 

(l + 7 2)(l +7 )2/3 
2 7 V3 (1 +7) 1/3 (7 _ 1) 

1 + 7 2 



ji2/3 



ji2/3 



(r, e,s) = T + 9T» + r Z7 "^ 1+7r : i7 - 1J r 2 / 3 + (s - r^)^±^-T^. 

(2.29) 

One might have noticed that in the scaling (I2.29P we extend x and y along 
the characteristic for r = 0, which is not exactly the characteristic for r ^ 0. 
However, this is not a problem, since the projection of (x(t, 6),y(r, 9)) along 



2/3 



7 



1/3 



the true characteristic on the line x + y 



(l+7) ; 



■T is (x(f,0),y(f,0)) with 



t = t + o(l). Then, the extension of Proposition 12.41 to the scaling (I2.29P is 
the following. 

Proposition 2.8. 

(a) Fix m 6 N, v G [0, 1), 77 G (0,1], and 7 G (0, 00) such that 
77> 7 (l + 7) -1 . Then, for any given real numbers t\ < t 2 < . . . < r m , 
61, . . . ,6 m and S\, . . . , s m , we have: 
(al) if 77 > 7(1 + 7 ) _1 , then 



YmvF ^]{L x {x{T k ,e k ),y{T k ,e k )) < £(r k ,e k , s k )}\ =P f]{A 2 (r k ) < s k } 
\fc=i / \fc=i 

(2.30) 

(a2) while ifrj — 7(1 + 7 ) _1 , then 

/ m \ / m 

Urn P ( pl^^Tfc,^),^,^)) < £(r fc ,4,s fc )} ) = P ( fl^BM^) < s fc } 



\fe=i 



\k=l 



(2.31) 
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(b) Fix m G N, i] G (0, 1] . Then, for any given 71 < 72 < • • • < 7m, such that 
rj < 71(1 + 7i)~\ and si, . . . , s m G R, ^ /ioWs 

f m ( f 2 1 \ 
MP H 7* W < ( ^ + ) # fc T + ^T 1 



\k=l 



rj 2 (1 — rj) 2 




<s k \\. (2.32) 



2.2.3 Two-sided LPP 



The main object of interest in this paper is last passage percolation with 
two-sided boundary conditions defined as follow. Given two paramaters 
7r, rj G (0, 1], the independent waiting times tUjj satisfy 

'Exp(l/7r) fori>l,j = 0, 



w 



■1 -j 



Exp(l/?7) 

Exp(l) 





if « = 0,j > 1, 

if i,j > 1, 

for i — 0,j = 0. 



(2.33) 



We denote by L 2 the last passage percolation time for the waiting times 
(I2.33p . where the subscript 2 stands for the two-sided. This corresponds 
with TASEP with two-sided Bernoulli initial conditions. The connection 
with two-sided Bernoulli initial condition for TASEP is obtained by setting 
rj = p_ and n = 1 — p + . 

The new phenomenon that occurs for two-sided LPP with respect to one- 
sided is the possible presence of Shockwaves in the corresponding TASEP 
picture. This occurs when tt + i] < 1, i.e., when characteristics meet. Indeed, 
the characteristic leaving from the axis (R+, 0) have slope (1 — ir) 2 ir~ 2 and, 
whenever rj + 71 < 1, they meet the characteristics leaving from the (0,R + ) 
axis, whose slope is r] 2 (l — t])~ 2 . The slope of the shockwave is determined 
by the Rankine-Hugoniot condition and it is given by the equation 

77(1 



7T 



y 



tt(1 -rj) 



x. 



(2.34) 



The limit shape is not anymore always as in f)2.32p but depends on the which 
side of the shockwave we focus on S(j) = lim^oo ^L2(T, r y 2 T) given by 

,,2 1 



7_ 
V 



+ 



77 



if7T<(l + 7)- 1 ,r / < 7 (l + 7)- 1 ,7 2 >^^ 



7 



1 - 7T 



+ Z> if 7 r<(l + 7)- 1 ^<7(l + 7)- 1 ,7 2 <^- 



TT 



(2.35) 
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When rj + n < 1, then the fluctuations are dominated by the boundary terms 
and live on a T 1//2 scale, while the bulk contribution to the fluctuations is only 
on a T 1//3 scale. Therefore, the limit process describing the fluctuations on 
each side of (not including) the shockwave is given by the Brownian motion 
obtained as the boundary contribution from the origin to the projections 
along the characteristics of the points we focus on. Thus, the two side of 
the shockwave will be independent. On the shockwave, there is a balance 
between the two boundary contributions: the last passage time for a point P 
on the shockwave is the maximum between the last passage time of the one- 
sided problem with w i = and the transposed one-sided problem u> j = 0. 
Since the fluctuations come only from the boundaries, the distribution of P 
will be the product of the distribution of the two one-sided problems, see 
Figure O 

This intuitive picture is confirmed by Theorem 12. 9[ which can be obtained 
from Proposition 12.81 without any additional hard analysis by using coupling 
arguments introduced in [TO] . 

Theorem 2.9. Consider the same scaling \2. 29\) as in Proposition ^. 81 
(a) Fix m G N, n, tc G (0, 1], 7 G (0, 00), and v G [0, 1). Then, for any choice 
of real numbers T\ < r 2 < . . . , r m , 9%, . . . , 9 m , and si, . . . , s m , we have: 
(al) Ifn > (1 + 7)^ and 77 > 7(1 + 7) _1 , then 

kpfp^(4 I 9 k ) lJ |(T tl fl i ))</(T kl 8M*)}) =p(f\{A 2 (r k )<s k }\ 

(2.36) 

(a2) If ir > (1 + 7) _1 and 77 = 7(1 + 7) _1 , then 



hm Pf f]{L 2 (x(T k ,6 k ) } y(T k ,e k )) < £(T k ,e k ,s k )}\ =P( f]{A B M^2(r k ) < s k } 

(2.37) 

(aS) Ifn = (l+ 7)- 1 and 77 = 7(1 + j)" 1 , then 

f]{L 2 (x(T k ,d k ),y(r k ,e k )) < £{r k ,e k ,s k )}) = p( f]{A stat (r k ) < s fe +r fe 2 } 

k=l ' ^fc=l 

(2.38) 

(b) Fix mi, m s , rrib G and set m = mi + m s + m& . zs i/ie number of 
points associated with characteristics from the left boundary of the LPP, mb 
to the bottom boundary, and m s on the shockwave (if exists). For ir, n G (0, 1] 
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and real numbers 9i,..., 9 m , choose j{ corresponding to each case. Then 
limpf f){L 2 (e k T,i*6 h T) < S{ lk )9 k T + suT 1 ' 2 } 

1 



k=l 
mi+ms 



7| ; 

t] 2 ( 1 — 77) ; 

1 7l 



7T- 



1-7T 



< s k 
< s k 



(2.39) 



where B and B' are two independent copies of Brownian Motion. 

Remark 2.10. The case tt — (1 + 7) -1 and rj > 7(1 + 7)" 1 can be recovered 
by Theorem 12.91 (a2) by the change of variable (x,y) — > (y, x) and 7 — > 7 -1 . 

Because of the nice correspondence with the TASEP, we focused here 
on LPP with two-sided "of width one", i.e., with modified weights only for 
i = and j = in f !2.33p . From a LPP point of view a natural extension 
is to consider weights different from 1 for a larger number of columns/rows. 
For example, the one-sided with boundary width equal to r was considered 
in [5]. Like for the Airy2 process j2], also one can describe joint distributions 
by PDE's PQ. For an extension to two-sided thick boundaries, see [T7]. The 
particular case of boundaries of sizes 1 and r corresponds, in terms of TASEP, 
to Bernoulli initial condition on Z_ with the first r particles having a different 
jump rate [16J. The coupling techniques would also work, except to the 
critical cases (like 77 + 7r = 1). 



2.3 Random sample covariance matrices 

We now give a random matrix interpretation of our result. For that 
purpose, we use a result proven in [22]. Consider an infinite array 
A(N) = [Aj] 

i>i,i<j<N where the Aij are independent complex Gaussian ran- 
dom variables with mean zero and variance l/(vrj + 7^). Define A(n,N) to 
be the matrix obtained by considering the first n rows from A. Define the 
NxN matrix M N {n) := A(n, N)*A(n, N) and denote by AAr )max (n) its largest 
eigenvalue. 

Let also L(n, N) be the last passage time from (1, 1) to (n, N) in the 
percolation model with independent exponential random variables with ex- 
pectations l/(7Tj + Tij). Then it is proved in [17J and [22] that the process 
n 1 — y Ajv,max(^) and the process of last passage times n h-> L(n,N) have the 
same distributions. Therefore, from Proposition 12.41 we deduce the following. 
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Theorem 2.11. Set 7Tj = 1 — 77, i > 2, 7i"i = 0, jtj — rj,j > 1. 

^ Fix m G N, 7 G (0, 00) and rj £ (0,1]. Recall x(r k ) , y(T k ), and £(r k , s k ) 

defined in IBM) . Set N = y(r k ) = jy^T . 

(al) If rj > j-p^ ; £aen /or real numbers T\ < t 2 < . . . < r m and si, . . . , s m it 
holds 



hm P ( f|{A JV>max (x(r fe )) < £(r k ,s k )} ] =P ( f^fo) < ) . (2.40) 



vi=l 



(ai?) If V = ^ en / or rea ^ numbers t\ < r 2 < . . . < r m and si, . . . , s m ^ 



holds 



lim P [ P|{A 



7V,max( 2<( T~k 



)) < ^(rjfc, s k 

)} = P ( n{^BM-^2(Ti) < S,} 



vfe=l 



vi=l 



(2.41) 

(b) Fix m £ N, rj £ (0, 1]. Se£ AT = |_^J , ^en /or any 71 < 72 < • • • < 7m 
snca taat r\ < 71(1 + 7i) _1 , and s±, . . . , s m £ M, it holds 

to p (q { WtJh - ("f + r~^) T + s * rl/2 }) 

= p(n{e 



vfc=l 



rj 2 (1 — n) 2 



<s k }\. (2.42) 



2.4 Limit processes: definitions 

Here we collect the definitions of the limit processes. 

Definition 2.12 (Airy 2 process, A2). The Airy-i process is defined in terms 
of finite dimensional distributions as 



L 2 ({Ti,...,r m }> 



(2.43) 



\k=l 



where Xs{ T k,x) — hx>s h ], an d is the extended Airy kernel: 



f dze^ T '-^ z Ai(s + z)Ai(s' + z), r > r', 
K A2 (r,s-y, S ') = { J % (2.44) 

dze (T ~ r)z Ai(s + *) Ai(s' + z) , t<t'. 
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The Airy 2 process was discovered in the PNG model [19]. It is a sta- 
tionary process with one-point distribution given by the GUE Tracy- Widom 
distribution F2 [55]. An integral representation of K^ 2 can be found in 
Proposition 2.3 of [36]; another form is in Definition 21 of [16] in the M = 

case. 

Definition 2.13. We denote by Abm^2 the transition process from Brownian 
Motion to Airt/2- It is defined in terms of finite dimensional distributions as 

P (C\{A B M->2(Tk) < s k } J = det(l - XsKA B u^Xs)L^{T lt ...,T m }xM), (2.45) 
where Xs{^k,x) = t[ x>Sk ], and K^ M ^ 2 is the rank-one perturbation Kj^ 2 : 

Ka^ 2 (r, s; t', s') = K M (r, s; r', s') + Ai(s) - J dze T ' z Ai(s' + z)\ 

(2.46) 

This transition process was derived in [33]. In Equation (3.6) of [33] the 
kernel is divided into two cases. However, using the identity (D.3) in 
namely, j R dye wy Ai(/3 + y) = e"^/ 3- ^ 1 ", we can rewrite 



^bm-j ( r > s > r '> s ') = k m ( r > s ; T \ s ') + Ai ( s ) / dze T ' z Ai ( s ' 



z) (2.47) 



for r' > 0. 

An integral representation of the kernel K^ BM ^ 2 can be found in 
Definition 21, in the M — 1 case. To see the Brownian Motion behavior, one 
has to take the r <C — 1 and replace s by s + r 2 . This shift is needed to take 
into account that the actual limit shape at the critical point changes from 
(I2.18P to (12.191) . So, for large — r, the approximation coming from (12.181) is 
not optimal anymore. Indeed, using (12.191) . s — r 2 would be replaced by s in 

The definition of the process for stationary TASEP, ^4 sta t, is quite intri- 
cate [6]. Its joint distributions is the r.h.s. of Equation (1.9) in [6]. 

Definition 2.14. The last process, B, is simply a standard one dimensional 
Brownian motion. Its finite dimensional distributions can be expressed in 
terms of a Fredholm determinant: let < T\ < • ■ ■ < T m , then 

P Hp{B(7fc) ^ Sk ^j = det ^ - X.#BX<Wn,wr™}xR), (2-48) 
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where Xs( T k, x ) = ^-[x>s k ], o,nd the kernel K B is given by 

AB(T ' S; T " ? > = ("27 J - V 2*(r - r-) < " P I" J ' 

(2.49) 



3 Proof of results 



3.1 Proof of Proposition 12.4 



Let n, N be positive integers. Consider the directed percolation model with 
one-sided boundary conditions f)2.17p . Let Li(n, N) be the last passage times 
from (0,0) to (n,N). The joint distribution of Li(n, N)jtl > can be ana- 
lyzed thanks to the so-called Schur process studied in |17) 4 I . In particular the 
joint distribution of the last passage times in the directed percolation model 
is given by: 

P ^p{Li(n fe , AT) < 5 fc }j = det(l-P s K N Ps) L * i{nu ... >nm}xm , (3.1) 
where Ps(k,x) := l[i>s fe ] and Kn is the correlation kernel given by 



K N {n i} x]nj,y) = -^„ ijnj (x, y) + K l N {n i} x;rij,y), 



K l N {ni,x;nj,y) = f dz f dw 



1 f I e wy ~ zx (z + 1- r])^- 1 z (w- r)) N 



v (2ni) 2 J c J c , w — z (w + 1 — r]) n i~ l w (z — n) N ' 
* ni , nj (x,y) = l [ni<n .]l [:c<2/] -L I dwe w{ y- x \w + 1 -7,)*-*. 

Z7T1 J c , 

(3.2) 

where C (resp. C) is a contour oriented anticlockwise and enclosing r) (resp. 
and r\ — 1). 



l The (N,p) in [T7] corresponds to (n — 1, N) in this paper. 
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3.1.1 The case where rj > 

Consider the asymptotics of the correlation kernel with the rescalin"' 5 

\ 2 / \ 2 

1 A T + _ ?^__ T 2/3 N 



l + 7 y (l + 7 )2/3 7 2/3 ' V X + 7 



(3.3) 



7 2/3 " 1 " l y/3 



We fix some s G K and assume that Sj > s for any i. To give the result we 
need a few definitions. Let us set 

7 (1 + 7) 1/3 p _ (l+7) 2/3 _ 1 + 7 (34) 

and Z(z) := exp (2rf/3 + 7}s$ + TFi(w c ) + s^T 1//3 pto c ) where 
2?iX\ _y_ 2 . / 1 27; 



(3.5) 

Proposition 3.1. Uniformly for s», s 3 - m a bounded interval, it holds 



Z(i) 

lim pT 1/3 ——K 1 N (n i ,x i ;n j ,x j ) 

N—>-oo ^\3) 



(3.6) 



/■oo 

= / e^"^ Ai(s< + rf + A) Ai(s,- + r 2 + A)dA + (^(T -1 / 3 ). 

Furthemore, for any k > 0, there exists a T large enough such that 

pT^^-K^x^x^ < Ce-« Si+s ^ (3.7) 

for all Si,Sj G M. and T > T . The constant C is uniform in T > T and 
Sj, Sj . 

Proof of Proposition \3.1\ The proof of Proposition 13.11 relies on a saddle 
point analysis of the correlation kernel (I3.2p with the rescaling (13. 3p . We 
first rewrite the singularity l/(w — z) in the kernel (13. 2p as 



1 



w — z 



e X( W -z) P T^ pT l/3 dX (3^ 







5 Comparing (|3.3p with (|2.20[) one sees two minor differences: (a) the integer parts are 
not explicitly written, but it is irrelevant for the large T asymptotics and (b) £ does not 
have the shift by — on s,. This is not a problem, because in our case the Tj are chosen 
in a bounded set. The scaling (|2.20p is reobtained in the end by replacing Si + t? by s,. 
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This allows us to rewrite as a product of two integral kernels: 

poo 

pT 1/3 K^(ni, Xi] Uj, xj) = / H(t u + A)G(r j , sj + \)d\, (3.9) 

Jo 

where 

-2n 



I |_r,^. T 2 / 3 

P7 



27T1 Jc (z-^y^i+y) 

Gfo,*) = e T^ xT ^ w l_ (w - yf^ 2 e s^ pwdw 



(w + 1 - 

(3.10) 

The contour C is like C but instead of anticlockwise it is clockwise oriented. 
Note that the two contours C and C still cannot cross each other. To perform 
a saddle point analysis of both G and H, we consider the first order leading 
terms in the exponential. Recall Fi from f)3.5p . Then for Tj in a uniformly (in 
T) bounded interval, Fi admits w c defined in (13. 4p as unique critical point. 
Furthermore 



i^> c ) = 2(l + 7 ) 2 ^T- 1 / 3 , F®(w e ) = -2^I±^1. (3.11) 

7P 7 

Now we briefly expose the ideas of the asymptotics since such arguments 
have already been developed many times (see e.g. [5] or at the beginning of 
the proof of Lemma 6.1 in [12] . where the steps are explained). From the 
assumption r\ > it follows w c > 0. Consider the contours 

d = {w c + te ni/3 ,teR}, C[ = {w c + te 27Ti/3 ,0<t<2}. (3.12) 
Then one has that 

^-ReFi ( Wc + = / + (7 - l)^-27t 2 t(t + l) 

dt l \ c 7 + 1 J (t 2 -t + l)(t 2 + 7 t + 7 2 ) tA l-t + t 2 

(3.13) 

The denominators in (13. 13j) is positive, being (product of) squared distances 
between the poles of the integrand and the integration path, e.g., t 2 — t + 1 = 
(7 + l) 2 |u> c + ^-e 271 " 1 / 3 + 1 — r]\ 2 . In particular ReFj decreases along C[. We 
now complete the contour in the upper half-plane as follows. Call wq the 
endpoint of C[ and set r := \w — 77 + 1|. Let < 9 < it be such that 
Wq = 1] — 1 + re l9 ° . Define 

C' 2 ■= {77- l + re ie A < 9 < vr}. (3.14) 
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0(7-1/3) 



Figure 3.1: The contours C and C are slightly deformed in the neighborhood 
of w c so that they don't touch each other. 



Then it is not hard to see that ReF; decreases along C 2 : 

^ReF t (v-l + re 10 ) = -r sin 6 (l + ^ - py^— ^) < ~cr sin 9, 

' (3-15) 

for some constant c > 0. Thus C — C[ U C' 2 U C[ U C2 is a steep descent path 
for Fi. 

For the z— contour, one has that 



ReF- u> -I e /3 = ■ — — - I-t-y - — T /3 

dt i V 7+1 ; (i + t + t 2 )(t 2 - 7 * + 7 2 ) i + * + * 2 ' 

(3.16) 

In particular -^ReFi(w c + ^-e 71 " 1 / 3 ) > as soon as t > j^jTiXT 1//3 . From the 
latter we deduce that the main contribution to the z— integral will come from 
a T -1 / 3 neighborhood of w c . The contour C\ is not a steep ascent contour 
for Fi but is enough for the purpose of evaluating the integral: it is a steep 
ascent path for the first order approximation of F iy that is forgetting for a 
while the 0(T _1 / 3 ) terms in Fj. 

In order to take care of the constraint on the contours C and C which can- 
not cross or touch each other, we now deform the w and z contours in a T -1 / 3 
neighborhood of w c so that the w- (resp. z-) contour lies to the left (resp. 
right) of w c (see Figure IO). As F[\w) = 0{T~ 1 ^), and F/ 3) (w) = 0(1) in 
such a neighborhood, the fact that the two contours are moved of 0(T -1 / 3 ) 
from the critical point has no impact on the asymptotics: this follows from 
a straightforward Taylor expansion of the exponential term. 

We then make the change of variables 



25 



Then it is not hard to see that for bounded s,, Sj it holds 

i -r /•ooe"' 1 / 3 

J_ e TF i (.w c )+. i( trV*v, cH , s .) = _L / e tV3s t t-r^ dt + 0( T -l/3) 

w c 2vri J^e-Tri/s 



-2T?/3-T iSi 



<D(T 



Wce -TF ] (w c )-s tP T 1 ^w, 



1/31 



(3.18) 

In the above we used Appendix A in [B] to derive Airy identities. In the case 
where s, > 0, one also gets the following exponential decay: let k' > be 
given. Then, as we can modify the contours C and C so that Ke(w — w c ) < 
— T i/ 3 while Re(z — w c ) > T i /3 one gets that 



w c e 

This ensures that 
lim pT 1,z ——K 1 N (n i ,Xi,n j ,x j 

N—too ^\J) 



-TFj{w c )-s iP T l ' s w c Qf Ti Si 



K Si 



< c- 



rpl/3 



< a 



K Si 



T l/3 ■ 



(3.19) 



poo 

/ e-^-^ Ai( Si +rf+A) Ai{ Sj +T*+X)d\ 
Jo 

(3.20) 

in the trace- norm class (one can choose k' = max{|rj|, i = 1, . . . , N} + k). □ 
We also need to consider the asymptotics of ///'' :i *•!>„.,. (./•,. .r ; * ' 



Proposition 3.2. For |s$ — sA in a bounded interval, it holds 



ffr 1/3 v n . tnj (x i ,Xj) 



i 



: exp 



Z(j) y^UTj - n) 
Furthermore for any k > 0, there exists a Tq large enough such that 



0(T- 1/3 ). 
(3.21) 



Z(i) 



< Qq~ K \ s i— s j\+(n-Tj)(si+Sj)/2 



(3.22) 



for all Si,Sj G E and T > T . The constant C is uniform in T > T and 



Si, Sj . 
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Proof of Proposition \3.2i The asymptotics of \& are again analyzed through 
a saddle point argument. Consider 

f(w) := X w - (P7)" 1 H™ + 1 - rj)- (3-23) 

Then 

M', i; ,J.r,. ) = _L£e 2 ^^^^ T2/3+ ^-^)- Tl/3 ^. (3.24) 

The critical point is again w c = r\ — and f"(w c ) = p 2 > 0. We choose the 
contour to be the circle centered at r\ — 1 and passing through u> c : this is a 
steep descent path. Making the change of variables w = w c + ^{ /3 one gets 
that 



27T Jr 



(3.25) 

Thus one obtains 

oT 1/3 ^ (x- x •le T(F<( ^ c)-Fj ' (Mc))+(si_Sj ' )Tl/3 ^ c+2(T ^ r ^ /3+riSi-riSj 

= /A = < exp + rf) + r lS% - r jSj ) + 0(T~^). 

^Ait{Tj - Ti) \ 4(r j -Tj) 3 J 

(3.26) 

The exponential decay for large Sj — Si is obtained as in [6], Lemma 21. For 
large positive Sj — Sj we can modify the contour so that it lies to the left of w c : 
Ke(w — w c ) < — n' /(pT 1 / 3 ) for any arbitrary k' > 0. For large negative Sj — Si, 
the contour is modified in the following way: it is again a circle centered at 
r) — 1 but passing through w c := w c + k' /(pT 1 / 3 ). It is a simple computation 
to check that / decreases along this modified contour and ensures that 

pT 1/3 ^ n . iinj (x i ,x j )e T(Fl(,i ' c) -^ (,i ' c))+(Sl -^ )Tl/3u;c < Ce~ K '\ Si - a ^. (3.27) 

The complete details to derive (13.221) from the above estimate is given in 
Lemma 21 in [6] (choosing again k' > max |tj| + k). □ 

Combining Proposition \'d.l\ Proposition 13.21 and the definitions of Sec- 
tion [53] yields part (al) of Proposition ET~4"1 
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3.1.2 The case where rt = 

' 1+7 

The rescaling is still given by ( 13. 3p . In this case, w c = 0. We recall that 
the contour C has to encircle the pole w — 0. Thus to get the exponential 
decay for large positive Sj, one needs to consider a different conjugation of 
the kernel. Indeed it is no longer possible to deform the contour C so that 
it lies to the left of w c . On the other hand, it is a well known fact that 
conjugation does not impact on the correlation functions of a determinantal 
random point process. 

Let then 5 > be given. Define 

Z{i, S) := exp Qr 3 + r lSi + TFi(w c ) + s lP T^ (w c + -^jA ) . (3.28) 
Proposition 3.3. Uniformly for s», Sj in a bounded interval, it holds 
Hm pT 1/3 Z Jf' 6 - K\ (m , Xi ; rij , x j ) 

= Ai( Sl + r?)e {s >- s i )5 (t 27 ' :i •'•'* - J™ Ai(r 2 + s 3 + x)e^ x dx^j 

POO 

+ e^-'t)* / e - x{n ~ T ^ Ai(si + r 2 + A) Aifo + r 2 + X)d\ + C(T- 1/3 ). 

(3.29) 

Furthemore, for any k > 0, there exists a Tq large enough such that 

P Tl/3 J^ K N(ni,Xi-,rij,Xj)\ < (3.30) 
for all Si,Sj G JR. and T > T . The constant C is uniform in T > T and 

Sj, Sj . 

Proof of Proposition VJ.'d Define then 

~ {pT^f r _ Tz _ 2r . xT 2 /3z (z + i- v )i^ +2T ^ _ s . T i /3pz , 

H(.Ti,Si) = . (be lz V '— 7 e s ^ pz dz, 

(3.31) 

We can now perform the saddle point analysis of the above kernels. We 
use the same contours as in the proof of Proposition 13.11 up to the following 
modifications: the contours C and C are deformed in a T _1//3 neighborhood 
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of w c = so that C encircles lying to the left of w c := w c + 5/ (pT 1//3 ) and C 
remains to the right of w c . Furthermore we can assume that the distance of 
these contours to w c is at least k' /T 1//3 . From the preceding, one easily gets 
for bounded Sj that 

i /•ooe"" 1 / 3 

Z(i, 5)H{n, Si ) = e^+K +T ^— / te t3/3 - s ^ t2 dt + 0(T~ 1 l' i ) 
= e^{-M'{^ + s l ) + T l ki{^ + s l )). 

(3.32) 

Similarly 

G{n,sj) =e - aiS -* l f/ 3 -* w J_ f°° e i e - t3/3+Slt+n ti dt + ( T -i/3^ (3>33) 

Z[i,d) 2vri yooe-ihri/s t 

where the contour passes to the right of 0. Using again Appendix A in [6] 

= e-^ f e -^/3^" - /°° Ai(rf + Si + x)e Ti *dx) + 0{T~ 1 ^ 



° (3.34) 
The exponential decay (as in ( I3.19P ) for large positive s» follows from the fact 
that the w— contour (resp. 2-contour) lies to the left (resp. right) of w c and 
with a distance at least /t'/T 1 / 3 . Again one shall choose k' > maxj |rj| + k. 

Finally, to derive the asymptotic correlation kernel in the case where Sj, Sj 
lie in a fixed bounded set, we use the simple algebra: 

e-^-)rfA™ / ds / dt(--V 



o 



(2vri) 2 yooe-iri/a V s / e s3 /3-^s-T^ 2 

ooe 27ri / 3 noce*' 1 / 3 ^ e * 3 /3-s r;i 2 ^ ' ' 



(27Ti) 2 iooe-avrt/s Jooe-.i/s S e s3 /3-s 3 s-T 3 s 2 ' 

yielding the asymptotic of given in Proposition 13.31 □ 

The asymptotic analysis of ^ for u> c = is almost unchanged from the 
last subsection. For bounded s« — Sj, one gets that 



( s >- s j) s ( ( Si - Si) 2 2 



(3.36) 
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To get the exponential decay (as in (13.221) ). one simply deforms the contour 
to the right or left of w c depending on the sign of Sj — Sj. 

Combining the above with Proposition 13.31 and definitions in Section 12.41 
yields part (a2) of Proposition I2.4[ by using the fact that 

det(l - X s ^iXs)L 2 ({T 1 ,...,T m }xii) = det(l - XsKiXa)v*({n,...,T m }xR), (3-37) 
where K\(s, x;t,y) = e b ^ x ~ y ^ Ki(s, x;t,y) for any b in a compact interval. 



3.1.3 The case where ri < -r— 

• 1+7 

Let m be a given integer. We now consider the asymptotic joint distribution 



(3.38) 



For k — 1 , . . . , m we set 
£fc = I — + 



r\ 1 — rj 



T + s k T 1/2 , n k = jlT, c k 



'll 



rf (1 — rj) 2 

(3.39) 

We recall that the correlation kernel Kn is defined in (13. 2p and let us set 



Z i) :-- 



for a given small 5 > 0. 



(l-r?r 

Proposition 3.4. For Sj — Sj in a bounded interval, it holds 
lim Vf^e^-^KNin^x^nj^j) 

N^oo Z(j) 

e (si-Sj)S / s 2 \ e {si-Sj)& 

3 eX P( -^5) - "[TKT,] 



7TC„ 



27T( C 2 - C? ) 



. exp 



2(c^ - cf ) 



Furthermore, for any k > 0, there exists a Tq large enough such that 



VT^e^~^ s K N ( ni , xf, n j} x 3 ) 



< Ce~ KlSi ~ Sj \ 



(3.40) 



(3.41) 



for all Si, s,; G K. and T > T . T/ie constant C is uniform in T > T and 



Si, Sj . 
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Proof of Proposition 3.4 Define 

, 2 



fj(w) := (^L + yz^J w + l) H w -V)~ Hw + 1 - rj). (3.42) 



KUn^x^n^xA = & dz & dw-^—e T ^-^ +Tl/2 ^ w ~ s ^. 



Then 

(2TTI) 2 J c J c , WW — Z 

(3.43) 

It is easy to check that the exponential term fj admits two critical points 
tx> c = and to" = ^^v^ < 0. The critical point w c = satisfies 

/ (0) = (1-77)2 — Tjj < 0. The steep descent (resp. ascent) path for fj should 
pass through w~ (resp. w c = 0). Nevertheless the contour C has to encircle 
the critical point w c . To deal with this difficulty, we separate the contribution 
of the pose at w = and the pole at w — 77 — 1. This will modify a little 
bit the saddle point analysis, which turns out to be similar to the analysis in 
Section 3 of [5]. 

Computing the residue at w — 0, one gets that 



K N (rii, Xi]7ij, Xj) = (p dz (b dw 



e WXj ZXi (z + 1 — 7]) T z (w — T])^ T 



' " ' (27ri) 2 Jc J c „ w - z (w + 1 - r/) T w ( z - n)^ T 

(z + l-r]) T (-v)^ T 



V "" -" ( ''" '' ?! (2ni)f c dZe ~ (1- V )T (z-r/)^' 

(3.44) 

where the contour C" does not encircle the pole w = 0. 
Let us first consider 

For ease we assume that 77 < 1/2 so that jj > 1, j = 1, . . . , m. Set the contour 
C 1 = { z = it, \t\ < 2}. Then 2Re(/<(W)) = 7 2 ln(t 2 + t? 2 ) - ln(t 2 + (1 - t?) 2 ) so 
that 

± Be( , , = . f _ 1 A = t (7, 2 -l)t 2 + 7Kl-^) 2 -^ 

dt KM )) \t 2 + ^ t 2 + (l-rj) 2 ) (t 2 + 77 2 )(£ 2 +(l-77) 2 ) 

(3.46) 

has the same sign as t. Thus Ci is a steep descent path for — fi with maxi- 
mum at z — 0. We complete Ci by the contour C 2 = {z = 77+ \/4 + 77 2 e ie , < 
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9 < 6> } where 9 is defined by rj + a/4 + r] 2 e W{> = 2it. Then setting 
z = 7] + y/A + ri 2 e w , 

^Mz) = "M*) (— + j— r— j » ( 3 - 47 ) 

a6* \ 77 1 — 77 |z + l — 

2 

and there exists a c > such that ^- + ^ — | 2+1 1 _ ?? | 2 > c along C 2 , so that 
also C2 is a steep descent path for — f\. 
Then, for bounded Sj we obtain 



"5^ 



lim y/f^rKj f (n i ,x i] n j ,x j ) = ^- [ e~* = -±^e 
n^oo Z{j) 2nd J R sjlnci 

(3.48) 

Next, we consider 



Kff(ni,Xi;nj,Xj) := <£ dz <f> dw 

{2m) z J c J c „ 



1 I j I j e wx i' ZXi (z + l-r]) T z (w-r])^ T 



y, w — z (w + 1 — rj) T w {z — r]) 7 i T 

(3.49) 

for some constant C > 0. We set C" = {\w~\e ld , 9 e [0, 27r]} so that w~ — r\ = 
— 7j(l — ?y) ensuring that 

^Re/ J (K|e ie )<-K|sin^S-^-^^) < -Csintf, (3.50) 

for < 9 < it. Thus Re/j achieves its maximum on C" at w = w~. Now a 
simple computation shows that fj(0) — fj(w~) = J w - f-(x)dx > 0. Thus, one 
has that 

Vfj^Kfrfaxtn^xj) <Ce~ cT , (3.51) 
for some constants C,c> 0. Thus, for bounded Sj we obtain 

lim v^f£^(n i} x i; n,-, Xj ) = 0. (3.52) 

To consider large positive s,, we consider the conjugated kernel 

^(n l ,x,;n i ,x i )e^-^ )<5 , (3.53) 

for some 5 > small. The contours Ci is modified so that it passes to the 
right of (5 + k) /T 1 / 2 . The w— contour passes to the left of (5 — k) jT 1 ! 2 . The 
exponential decay for large positive s$ follows. 
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The analysis of ^f ni ,n is similar to those of the preceding sections. The 
exponential term to be considered is 

g(w) := ( 7 , 2 - t|) (jj + Hw ~ v)) ■ (3-54) 

2 2 

g has a single critical point w c = with g"(0) = — l —^ J - = c 2 — cf. Consider 
the contour w = r/(l + e ie ). Then the leading (i.e., non exponentially neg- 
ligible) contribution to *& ni ,nj comes from a neighborhood of width T -1 / 2 of 
w — 0. Thus for bounded |sj — s»| we deduce that 

lim > /r|®-* n . ft .(a; i> a; i ) = — / dt e -( c ?- c ?)* 2 /2+i^- S ;) 
r-^oo Z(j) J ' * ; 2vry M 

1 ( {si-Sj) 2 \ ( 3 - 55 ) 



eX P \ o / 2 2 x 

27r( S 2 - C ?) V 2(c 2 -c 2 

For large |sj — Sj|, we consider the conjugated kernel 

ir JV (n i ,Xi;n i ,x i )e (Si_ ^ ) ' 5 . (3.56) 

Depending on the sign of s, — Sj we modify the contour to be the circle of 
ray 5 ± 2k' / a/T. This ensures the exponential decay for large | Sj — Sj\ . □ 

The above Proposition 13.41 has the required asymptotic results needed to 
conclude part (b) of Proposition 12.41 



3.2 Coupling Lemmas 

The following technical lemmas provide a basis for the coupling arguments 
necessary in our proof of Theorem 12.91 They provide generalizations of 
Lemma 4.1 and Lemma 4.2 of [10J from one-point functions to n-point func- 
tions. The proofs, however, are almost identical. 

For the purpose of the lemmas let X n and X n take values in M fc . We 
say that X n > X n if, with probability one, every coordinate of X n is greater 
than or equal to the corresponding coordinate of X n . We say that X n =3- F, 
where F is a distribution function on R fc if for all e > and (si, . . . , s&) a 
continuity point of F, there exists an iV = N(e, si, . . . , such that for all 
n > N, 

|P {ntAK < Si}) - F( Sl , . . . , s k )\ < e. (3.57) 

Finally we say that X n — X n converges in probability to zero if for all e > 0, 
lim^oo P(||X„ — X n ||oo > e) — > (the infinity norm is just that max over all 
finitely many coordinates: ||^||oo = niaxi<j<fc \X l \). 

In the three lemmas below we assume all random variables are M. k valued. 
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Lemma 3.5. If X n > X n and X n =^> D as well as X n =>■ D, then X n — X n 
converges to zero in probability. Conversely, if X n > X n , X n =>■ D and 
X n — X n converges to zero in probability then X n =>• D as well. 

For vectors X and Y in R fc we define Z = max(X, Y) to be the coordinate- 
wise maximum (i.e., Z % = max(X J , Y 1 ) for i — 1, . . . , k). 

Lemma 3.6. Assume X n > X n and X n =>■ D\ as well as X n =>■ D\; and 
similarly Y n > Y n and Y n =>- D 2 as well as Y n =>■ D 2 . Let Z n = max(X n , Y n ) 
and Z n = max(X n , Y n ). Then if Z n D^, we also have Z n D%. 

Lemma 3.7. Assume X n > X n and X n Di as well as X n D\; and 
similarly Y n > Y n and Y n =^ D 2 as well as Y n =>■ D 2 . Then if (X n , Y n ) =>■ F 
(a 2k- dimensional distribution function) so does (X n ,Y n ) F. More gen- 
erally this also applies to m sequences of random variables under the same 
hypotheses on each sequence and on their m-point joint distribution function 
limit. 

Proof of Lemma Iff. 51 This proof is a straight forward generalization of the 
proof of Lemma 4.1 of [10] and hence we will not reproduce it. All inequalities 
in the original proof should now be considered as holding true coordinate- 
wise, all absolute values should be replaced by norms (on M. k ). The e sized 
blocks used should be replaced by k dimensional e boxes, and all intervals 
should interpreted as boxes in R fc . Other than these changes, the proof goes 
through word for word. □ 

Proof of Lemma Iff. 61 Again this proof is word for word the same as Lemma 
4.2 of [10], with the modified interpretations of notation noted above. □ 

Proof of Lemma |ff. 7[ Lemma 13.51 shows that X n — X n and likewise Y n — Y n 
converges in probability to zero. This implies that for all e > 0, using the 
triangle inequality and the union bound, 

P(\\(X n ,Y n )-(X n ,Y n )\\ 00 >e) < P(\\X n -X n \\ 00 >e) + P(\\Y n -Y n \\ 00 >e) 1 

(3.58) 

which, by Lemma 13751 goes to zero as n — > oo. This immediately implies that 
the joint (X n , Y n ) converge to the same distribution as (X n , Y n ). Lemma 
is, in fact a corollary of this result. The generalization follows by the exact 
same argument as above. □ 
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3.3 Proof of Theorem MM 

In Section 13.11 we proved Proposition 12.41 directly from asymptotic analysis 
of the Schur Process. From that theorem we will, using the three lemmas 
above and the slow decorrelation result of Proposition 12. 5[ provide proofs of 
Proposition 12.81 and Theorem 12.91 

Proof of Proposition \2.8l The proof is based on Proposition ^. 4l together with 
slow-decorrelation phenomenon (Proposition 12.51 see also Remark 12. 7p . 

Consider first cases (al) and (a2). Denote by f k the number such that 
(x(h),y(h)) (defined in ^M) and (x(r k ,9 k ),y(r k ,9 k )) (defined in (gM) 
belongs to the same characteristic line. Moreover, notice that the projection 

2 

along the characteristic direction for r = to the line y = jpp^T is obtained 
by choosing 

in (I2.29p . However, the slope of the characteristic line passing by 
(x(r k , 9 k ),y(r k , 9 k )) differs from the slope of the characteristic line for r k = 
by just 0(T -1 / 3 ). Therefore, asv<l,f k =r k + O^T"' 1 ) -> r fc as T ->■ oo. 
Also, due to slow decorrelations, the fluctuation of Li(x(r k , 9 k ),y(r k , 9 k )) dif- 
fers from the fluctuation of Li(x(f k ) , y(f k )) by o(T 1 / 3 ), whose differs from the 
fluctuations of Li(x(T k ),y(T k )) again by o(T 1 / 3 ). Thus Proposition (al) 
and (a2) follows. 

The case (b) is even simpler. In that case, the point (9 k T, ^l9 k T) is on 
the same characteristic line as (T, j k T), at a distance 0{T). Therefore by 
Proposition 12.51 (b) the fluctuations of Li(9kT,jl9 k T) and Li(T, J^T) differs 
only by o(T l l 2 ), from which Proposition 12.81 (b) follows. □ 

Proof of Theorem \2.!A We follow the method of [TD] and define two coupled 
random vectors X and Y. X is the vector of last passage times from (0, 0) 
to (x(Ti,9i),y(Ti,9i)), 1 < i < m, with last passage paths forced to take a 
first step to the right, and Y is the vector of last passage times with paths 
forced to take a first step up. Therefore, their coordinate-wise maximum 
Z = max(X, Y) is the last passage times without any restrictions on the first 
step (i.e., Z % = L2(x(T i ,9 i ),y(T i ,9 i )). A key observation is that X and Y 
are both marginally distributed as the last passage times for last passage 
percolation models with only one-sided boundary conditions (as opposed to 
the two-sided conditions we must consider for Z). In the case of Y the one- 
sided boundary waiting time is exponential of mean 1 / 77 and in the case of X 
the boundary waiting time is exponential of mean 1/tt. But the coordinates 
must be flipped so as to conform to our definition of last passage percolation 
with one-sided boundary conditions (the boundary condition should appear 
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on the left boundary, not the bottom). Depending on the regime of fluctua- 
tions we will be able to compare the random vectors to related but simplified 
vectors A and Y which have the same asymptotic limiting distribution but 
are strictly less than A and Y. Then, using the coupling lemmas we will be 
able to show that Z and Z = max(A, Y) have the same distribution limits as 
T goes to infinity. However, as X and Y are simpler than X and Y, we will 
be able to identify Z exactly and hence determine the asymptotic multipoint 
distribution of Z, completing our proof. 

As slightly different coupling arguments are necessary for each part of the 
theorem we will split the proof up according to the four cases of the theorem. 



3.3.1 Proof of Theorem [Ml part (al) 

This case corresponds to both X and Y being last passage time vectors from 
last passage percolation models with one-sided boundary conditions of small 
enough mean so as to behave asymptotically the same as the corresponding 
models without boundary conditions. With this in mind we define X and 
Y, random vectors which are coupled to X and Y in terms of the underlying 
random last passage waiting times. Let X = X and let Y be the vector of last 
passage times defined by Y l = Li(x(Ti,9i),y(ri,9i)). The new last passage 
time Li is the last passage time in a coupled model where the boundary 
waiting times (which are exponential with mean 1/rj) are multiplied by rj 
(hence making them distributed as exponentials of mean 1). The key is 
that this random last passage time is coupled to the last passage time L\) 
since that they are based off of the same random waiting times. Additional, 
because of r\ < 1 it holds L\ < L\. Therefore Y < Y where the inequality 
is in terms of each coordinate separately. More trivially we also have that 
A < A. 

In order to apply our coupling lemmas we must center and rescale X, Y, X 
and Y so that our new X = (A 1 , . . . , X m ) equals the vector with coordinates 

A*-l(r^,0) 
( l + 7) 2/3 7 -i/3 r i/3 ' !<*<™> (3-60) 

with £{t, 8, 0) given in (I2.29p . and likewise for the other variables. Under 
this centering and rescaling A < A and Proposition 12.81 shows that both 
A and A converge in joint-distribution to the same Ai process as T — > oo. 
Likewise Y < Y by construction and Y and Y converge in joint-distribution 
to the same A<i process as well. Moreover, Z = max(A, Y) is (except for a 
single waiting time of zero at the origin, which is asymptotically irrelevant) 
the last passage time vector for a one-sided last passage percolation model 
with boundary waiting times with mean l/rj. Proposition 12.81 shows that Z 



36 



converges in joint-distribution to the A2 process. Therefore, using Lemma I3TB1 
it follows that Z = ma.x(X, Y) also converges to the A2 process, which is 
exactly what we needed to prove. 

3.3.2 Proof of Theorem [2J2 part (a2) 

We are in the case of r\ — 7(1 + 7) _1 and n > (1 + 7)" 1 . As such, we can 
apply the exact same argument as in the proof of part (al) above. The only 
difference is that the Z process will now, as determined by Proposition 12.81 
converge to the *4bm->2 process. Therefore Z will also converge in finite- 
distribution to the Abm^2 process, which is, again, what we desired to show. 

3.3.3 Proof of Theorem [2J2 part (a3) (see [6]) 

This proof is the subject of the recent paper [6]. The coupling techniques 
employed for all of the other proof do not apply here. The heuristic explana- 
tion is that the last passage time comes from the competition of two sets of 
paths each of which goes along the boundary for distance of order T 2//3 and 
then enters the bulk. Because the range of the transversal fluctuations of a 
last passage path are of that order T 2 / 3 , these sets of paths have non-trivial 
correlation, which is evident in that fact that they yield a different process, 
the A^t process. 

3.3.4 Proof of Theorem [Ml part (b) 

As before we write the last passage random variable L2(6iT, / yf8iT) as 
max(X l ,Y l ) where X % and Y l are coupled last passage times, restricted to 
paths which step first right or up, respectively. We now couple X 1 with X 1 
which is the last passage time when forced to stay along the bottom edge for 
a specific deterministic fraction of the path, and then depart into the bulk. 
Specifically we define X % to be the max of passage times over all paths which 
go distance 



and then take a step up. Likewise we define Y l to be the max of the passage 
times over all paths which go distance 



and then take a step right. It is clear that X 1 < X % and that Y % < Y l . What 
is not obvious is the choice of distances. In short, this is given by the solution 




(3.61) 




(3.62) 
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to an optimization problem at the level of the law of large numbers (see [3] 
for an explanation of this heuristic). 

Define the following events: for i G {1, . . . , m/ + m s }, set 



T] 1 — T] 

while for i G {1 + mi + m s , . . . , m}, set 



l T + s t T 1 ' 2 



(3.63) 



E; 



L^T^T) < ( - + 



7f 



7T 1 — 7T 



6;T + s,T 1/2 



(3.64) 



Notice that for i G {m^ + 1, . . . ,mi + m s } both definitions are identical. 
Likewise, define Ei except in place of L 2 (9 i T,'y 2 9iT) use Z % = max(X*,F*). 
Let us denote L 2 (9 i T,j 2 6 i T) = max(X i ,Y i ) as Z\ It is clear that Z* < Z\ 
We claim that 



T ih„p(n^) = itap(n^) 

v fc=i 7 v fc=i 7 



(3.65) 



For % G {1, . . . , mi} center and scale X 1 , X 1 , Y l , Y l , Z l and Z l by applying 



x 



x i — y 



v ^ i-v) 



1*1/2 



(3.66) 



It follows from the one-point fluctuation result of [10] that the centered and 
scaled Z % and Z l both converge in distribution to the same Gaussian random 
variable with variance 



0; 



7/ 



if (1 — rj) 2 

For i G {mi + m s + 1, . . . , m} we center and scale with 



(3.67) 



x 



x i — y 



A + _Zl_ \Q T 

7T 1 — 7T / ' 

TV2 



(3.68) 



Then the centered and scaled Z % and both converge in distribution to the 
same Gaussian random variable with variance 



6; 



1 



7T 



Of 



1-7T 



(3.69) 



For i G {m/ + 1, . . . , mi + m s } we center and scale as (I3.66p . Then the 
centered and scaled Z l and Z % both converge in distribution to the maximum 
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of two Gaussian random variables with variances given by equations ( I3.67P 



Since for every % £ 1, . . . , m, the centered and scaled Z l and Z % converge 
to the same distributions, and since Z % < Z l , Lemma 13.71 implies that the 
asymptotic joint distribution of the Z l , and of the Z % converge to the same 
distribution. This proves the claim given in equation (I3.65p . 

Therefore it remains to show that the joint distribution of the Z % behaves 
as desired. The fluctuations given by Z % are a combination of the fluctuations 
from the boundary waiting times and from the bulk waiting times. However, 
since we scaled by T 1 / 2 and since the boundary and bulk fluctuations are 
independent, the bulk fluctuations have a prefactor of T -1 / 6 and hence (by, 
for instance applying the Converging Together Lemma on page 89 of [24J) 
only the boundary fluctuations contribute asymptotically. The covariance 
of these fluctuations depends on portion of the boundary which the X 1 and 
Y % depend upon. We can encode this covariance structure in terms of two 
independent Brownian motions (one for the left boundary and one for the 
bottom boundary). For i £ {1, . . . , m{\ all of the fluctuations come from the 
left boundary. For i £ {mi + m s + 1, . . . , m} all of the fluctuations come from 
the bottom boundary. For i £ {mi + 1, . . . , mi +m s } fluctuations come from 
the maximum of the left and bottom boundary Brownian motions. Writing 
down this joint distribution leads exactly to (I2.39p . □ 

Remark 3.8. It is worth noting that in the proof of part (b) above we 
did not, in fact, appeal to the analogous one-sided last passage percolation 
result of Proposition 12.81 This is because we needed to establish the product 
structure and hence reduce everything to just processes along the boundary. 
As such the Brownian motion results of Proposition 12.41 and 12.81 may, in fact, 
be proved directly in this manner (as they are corollaries of this result) and 
do not require the asymptotic analysis of the Schur process. 

3.4 Proof of the TASEP height function theorem 

Proof of Theorem \2.1[ The connection between two-sided directed percola- 
tion and TASEP has been discussed in Section 12.2. II Consider first the cases 
(al)-(a3). From (I2.16P we have 



and (|gm 




k=l 




k=l 
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with (difference of order 1, due to the integer parts which are not explicitly 
written since they are irrelevant in the asymptotics) 

x k = ~pr(r fc , 6 k ) + H (r fe , 9 k , s k )), ^ 
y k = ±(H(r k ,O k ,s k ) - X(r k ,O k )), 

where X and H are defined in (I2.9p . Explicitly, by setting 7 := (1— £)/(l+£), 
i.e., £ = (1 — 7)/(l + 7), we have 

** = 7^^7(T + fc T") + r fc / / 7 - T 2 / 3 + (r 2 - s fc W 7 ^ T 1 / 3 , 



2/fe 



(1 + 7 ) 2V ,K (l + 7 ) 5 / 3 Vfe K7 (l + 7 ) 4/3 ' 

" " iT + k T v ) - r-^—T 2 / 3 + (r 2 - s^—^—T 1 '*- 



(1 + 7 ) 2V 7 (l + 7 ) 5/3 V " K7 (l + 7) 4/3 ' 

(3.72) 

Once the problem is rewritten in terms of directed percolation, the the- 
orem is proven using Theorem 12.91 and the slow decorrelation (see Proposi- 
tion !2.5p . i.e., we use a similar strategy of the proof of Proposition 12 .81 starting 
from Proposition 12.41 

For the above given (x k ,y k ), the limit shape ( I2.18P gives us 



(I 1 ^2/3 

x k (l + ^y~k~/x~ k ) 2 = T + k T v - s k K ^ T 1 / 3 + 0(1). (3.73) 



Therefore, 



p(f]{L 2 (x k ,y k ) <T + O k T»}\ 

f| {L 2 (x k , y k ) < x k {\ + y/y k /x k ) 2 + s k - ^— T 1 ' 3 + 0(1)} J . 

fc=i ^ ' 

(3.74) 

The fluctuations (with respect to the limit shape behavior) are, by the slow 
decorrelation theorem, the same as the fluctuations of the projection along 
the characteristic line on the line x + y — t-^yi (T + 8 k T u ). Exactly as in 
the proof of Proposition 12.8} we can use an approximate characteristic line, 
namely the characteristic line for 77 = s k = 0. We look for f k such that 

x k = x(f k , k ) + r(l + 7)~ 2 
Vk = y(h,o k ) + r7 2 (l + 7) 



2m ..,-2 (3-75) 



with x(T,0),y(r,0) as defined in (I2.29p . If f k — > r k as T — > 00, then 
the theorem is proven. This is the case, algebraic computations lead to 
f k = r k + 0((s k - r^T- 1 / 3 ) as desired. 
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Consider now the case (b). From (I2.16P we have 

f]{he k T^k0 k T) > h^k)e k T-2s k T l ' 2 }\ = p( f]{L 2 (x k , y k ) < 6 k T} 

k=l ' ^fe=l 

(3.76) 

with 

x k = \{£ k 6 k T + h ma (Z k )6 k T - 2s k T^ 2 ), 
y k = \{h^ k )6 k T - i k 9 k T - 2s k T 1 / 2 ). 

Let us focus on the case < 1 — (p_ + p + ) (remind r\ = p_ and n = p + ); we 
have 

x k = (1 - p_)(p_ + &)0 fe T + s k T 1 ' 2 , 
y k = p_(l - p_ - ^ k )6 k T + SkT 1 ' 2 . 

Let and 7^ such that x k = 6 k T and y k = O^lT. Then, we get 

T 1 / 2 

9 k T = S(%)6 k T + s k — -, (3.79) 

P-(l-P-) 

and 

4 = (1 - p_)(p_ + 6)4 + ^(T- 1 ^). (3.8O) 

Then, we can apply directly (12.391) with 9 k replaced by 9 k , ^ k by 7^ and s k 
replaced by s k / (p-(l — p~)) (similarly for the case > 1 — (p_ + P+)) to get 

limP( n{L 2 (x fc ,y fc ) < 9 k T} 

T— too \ 1 < 



fc=l 



fc=l 



p ( n {* fa 



p-(i-p-) ; - p-(i-p-) 

£ fc + 2p + - 1 \ s k 



< 



p+0--p+) J p+i 1 - p+) 



\ k =mi+l 

Finally, one uses the scaling of Brownian Motion to rewrite it as in (I2.14p . □ 
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